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Abstract: 

It has been suggested recently that an arbitrary induced ^-vacuum state {6 md ) could 
be created in heavy ion collisions. If such a state can be created, it would decay by 
various mechanisms to the Qf und = state which is the true ground state of our world. 
In the following we will discuss the possibility of studying this unusual state through the 
emission of pions, ^-mesons, and //-mesons. We will also present the spectrum of the 
produced particles in this non-zero 9 background. We use the instantaneous perturbation 
theory for our estimates. 



1 Introduction 



The long awaited Relativistic Heavy Ion Collider at Brookhaven (RHIC) could yield many 
fascinating results concerning the accepted theory of the strong force, QCD. Although the 
main goal of this experiment is to study the quark-gluon plasma, there exist many other 
possibilities. One of the most sought after results is the nature of the chiral symmetry 
restoration and/or deconfining phase transition. Another possible outcome is the so called 
disoriented chiral condensate (DCC). In addition, it was recently suggested U that it may 
be possible to create an induced 9 7^ 0-state in heavy ion colliders. If this is possible, the 
consequences could include the observation of a new state of matter. 

Indeed, as is known, in the infinite volume limit and in thermal equilibrium the 9- 
vacuum state is the absolutely stable ground state of a new world with new physics 
quite different from ours. In particular, P and CP symmetries are strongly violated 
in this world. In spite of the fact that the ^-vacuum state has a higher energy (see 
below) this state is stable due to the super selection rule: There does not exist a gauge 
invariant observable A in QCD, which would communicate between two different worlds, 
(9'\A\9) ~ 5(9 — 9'). Therefore, there are no transitions between these worlds |J. 

Of course, we do not expect that such a stable #-state can be produced in heavy 
ion collisions. However, we do expect that locally, for a short period of time due to the 
non-equilibrium dynamics after the QCD phase transition, a large domain with wrong 9 
direction may be formed. This provides us with a unique opportunity to study a new 
state of matter. The idea is very similar to the old idea of the creation of the Disoriented 
Chiral Condensate (DCC) in heavy ion collisions || -||. DCC refers to regions of space 
(interior) in which the chiral condensate points in a different direction from that of the 
ground state (exterior), and separated from the latter by a hot shell of debris. Our starting 
point is the conjecture that a classically large domain with a nonzero U(1)a chiral phase 
may be formed in a heavy ion collisions, i.e. we assume that the expectation value for the 
chiral condensate in a sufficiently large domain has, in general, a non-zero U(l) phase: 
(^l^r) ~ e l ^ a '" 3let \(^/^/)\. This phase is identified with 9 md . Such an identification is a 
direct consequence of the transformation properties of the fundamental QCD Lagrangian 
under U(1)a rotations when the chiral singlet phase can be rotated away at the cost of 
the introduction of the induced 9, 9 md = —Nf<p smglet . 

We should remind the reader at this point that in QCD, the fundamental 9 parameter 
enters in the combination 9Q = 9^G a G a . This gives us the first contribution to the 
physically observable parameter 9. The second contribution is related to the phase which 
is introduced upon diagonalization of the quark mass matrix, which can be in principle 
arbitrary. Current experimental results indicate that these two numbers cancel with 
precision better than 10~ 9 . The problem of why 9 ~ in our world is known as the strong 
CP-problem, and we do not address this problem in the present paper. Our point is the 
following: by having a non-zero 17(1) phase {V L V R ) ~ e**™" 1 " | | in a sufficiently 
large region, we essentially break this precise cancellation between two contributions in 
a macroscopically large domainQ. Exactly this non-cancellation will mimic a non-zero 
^-parameter. 

■""We use the term macroscopic in this paper to refer to scales which are larger than the fundamental 
scale, A QCD 
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The production of non-trivial # md -vacua would occur in much the same way as dis- 
cussed previously for DCC[||-||. The new element is that, in addition to chiral fields 
differing from their true vacuum values, the induced ^-parameter, which is zero in the 
real world, becomes effectively nonvanishing in the macroscopically large domain as ex- 
plained above. Therefore, although the fundamental ^-parameter remains zero, we would 
expect that an induced ^-parameter can be non-zero in a macroscopically large domain, 
and it will mimic the physics of the world with Qf und ^ 0. A different name for the 
same object would be a " zero (spatially independent) mode of the rf field" which ex- 
actly corresponds to the flavor-singlet spatially independent part of the chiral condensate 
phase. However, we prefer to use the term 9 md because the symbol rf is usually associated 
with the free asymptotic state of the heavy 7/-meson (excitation) and not with a classical 
constant field (condensate) in a large domain where r/'-mesons interact very strongly and 
change their properties in this background. This other terminology could be also mislead- 
ing because of this reason. In addition to the above mentioned difference, the disoriented 
chiral condensate involves the amplification of the low momentum modes of the charged 
pions, while the formation the 0-state involves the amplification of the low momentum 
modes of only the neutral particles, including the rf singlet . 

One might ask, "How could such a md -state be detected if it is created in heavy 
ion collisions for a very short period of time?". One obvious possibility is that the Q ind - 
state could be observed through Goldstone bosons with specific CP-odd correlations 0. 
However, a more promising signature is related to the low momentum modes of the 
chiral fields 0. Indeed, as was shown in [I], the creation of the mrf -state involves the 
enhancement of the low frequency modes of the neutral chiral fields. Therefore, it is 
natural to expect that particles made of these chiral fields are to be produced when the 
9 m -state blows apart. To be more precise, in the limit m u = = m s <C Aqcd the 
only relevant physical degree of freedom is a singlet phase of the chiral condensate which 
transforms, at the very end of the transition exclusively, to the ?/-mesons when the Q md - 
state falls apart. Therefore, in this limit we would expect an excess of r/-mesons to be 
produced. However, due to the mixing and difference in masses, m u 7^ irtd 7^ m s , it is 
expected that all neutral modes are to be produced. The main goal of the present work 
is to study the spectrum and the angular distribution of particles which will be produced 
when the md -state falls apart. 

However, before we discuss some numerical estimates, we would like to present some 
estimates based on simple energetic considerations. As we discuss below (|]), the vacuum 
energy density in the 0-state is greater than in the 9 = vacuum state by the amount: 



where for simplicity we use Nf = 2 with m u = rrid = m. When the 9 -state blows apart, 
the energy associated with this background will be released as Goldstone bosons carrying 
the total energy: 



where V = L 3 is 3d volume of the # md -background measured in Fermi. Therefore, given 



AE = E vac (9) - E vac {9 = 0) = 2m|(gg)|(l - I cos-l), 
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the cm energy y/s = 40 TeV, only a small fraction: 



AEV 20( W ^ F ) GeV _ J y 



p 10 



y/s 40 TeV \(10 fmf J 

of the total collision energy will be released through the decay of the induced # mci -state. 
Two of the most profound features of the Goldstone bosons produced in this decay are 
the following: 

1. Due to the fact that the 9 md formation is caused by amplification of low momentum 
modes, we expect the spectrum of the Goldstone bosons from this decay to be strongly 
enhanced at low \k\ ~ L -1 ; 

2. Due to the nonzero value for the topological density (6\ffcGG\B) = - dVv gg {9) = 
— m\ (0|gg|0) | sin |, see (|6|), one could expect that some P, CP odd correlations would 
appear in this background. 

Considering this, we should look for enhanced production of particles such as the 
71,1], and t]' with low momentum, on the (10 — 100) MeV scale (depending on the size 
of the domain L). The difference between the vacuum state and the DCC which also 
involves enhancement of low momentum modes is that the ^-vacuum state also produces 
//-mesons. These should decay eventually by various processes to photons and dileptons, 
and if a large number of these low momentum particles are detected this would be a 
definite signal of the creation of a macroscopically large domain with 8 md ^ 0. One 
could further speculate that such an enhancement could also account for the unexplained 
large number of low momentum dileptons seen at CERN || as suggested in[[|. However, 
detailed calculations are needed to give this speculative conjecture further support. 

The presentation is organized as follows. In Sect. 2 we briefly describe the effective 



Lagrangian constructed for QCD as given in [Tjj, [TT[| . In Sect. 3 we demonstrate that it 
is theoretically possible to create 8 md ^ 0-states in heavy ion collisions. The spectrum 
of the Goldstone bosons produced in the nonzero 9 background is calculated in Sect. 4 
for different geometries of the # md -region. Finally, we end with concluding remarks and 
future considerations in Sect 5. 



2 Low Energy Effective QCD Lagrangian 



The starting point of our analysis is the low energy effective Lagrangian which reproduces 
the anomalous conformal and chiral Ward Identities. The corresponding construction in 
the large N c limit has been known for a long time[|12|, |i~3f . 

for finite N, 



construction of ref. 12, 13 



10, 11 



The generalization of the 
and we shall use formulae 



from the papers |T0, ITT 



was given m 

However, we should remark at the very beginning that all local 
properties of the effective Lagrangians for finite and infinite N c are very much the same. 
Small quantitative differences in local physics between the description of |12|, |Tj| on the 

|TT|j on the other hand, do not alter the qualitative 



one hand and the description of []T0|, 
results which follow. 

The effective Lagrangian describes the light matter fields of QCD which consists of 
an octet of pseudo-Goldstone bosons (7r's, K's, and the if), and the r( singlet field. We 
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parameterize these fields by a unitary matrix Uy corresponding to the 75 phases of the 
chiral condensate ((^r^l) = — \{^R^L)\Uij) in the following way: 



U = Uq exp 



r ir a \ a 2 rf 
IV 2— ■ Vi- 



fn */Nffl 



f Jr l 



where Uq solves the minimization equation for the effective potential, 7r a represents the 
pseudoscalar octet, and Nj is the number of flavors. The effective potential derived in 
|r0| , [TT] ] takes the following form: 

V(U,0) = -Ecos[— (9-ilogDetU)] - ]-Tr(MU + AfilP), (2) 

iV G A 

where M is the diagonal quark mass matrix defined with their condensates: M = 
-diag(m i (W i ^ i )) and E = (ba s / (32tc)G 2 ) is the vacuum energy with b = HJV c /3-2iV//3. 
Expanding the cosine (this corresponds to the expansion in 1/N C ) we recover the result of 
|12|| at lowest order in 1/N C together with the constant term E required by the conformal 



anomaly: 



V(U, 9, N c -> 00) = -E - ~~ /YM (9 - i log Bet U) 2 - -Tr (MU + M + U + ) +... . (3) 

Here we used the fact that in the large N c limit E(-^-) 2 = —{u 2 )ym where {u 2 )ym < 
is the topological susceptibility in pure YM theory. Corrections in 1/N C stemming from 
Eq.(^) constitute a new result of ref.[10|. In order to demonstrate that the effective 
potential shown above is of the correct form, we can verify that the following three features 
are characteristic of Eq.(Q): 

1. Eq.dH) correctly reproduces the Witten-Di Vecchia-Veneziano effective chiral La- 
grangian 0] in the large N c limit; 

2. It reproduces the anomalous conformal and chiral Ward identities of QCD; 

3. It reproduces the known dependence in 9 (ie. 2n periodicity of observables) [ |T2|j . 

To study the vacuum properties of the 9 vacua it is convenient to choose a diagonal 
basis to parameterize the fields U as U — diag(expi<j) q ), q = u,d,s such that the potential 
(0) takes the form: 

V(M<) = ~E cos( 9 ~^ 0i ) -J2 M i cos «fe ( 4 ) 

where Mj are the diagonal entries of the quark mass matrix which was introduced above. 
Notice that in Eq.(j4|) the 9 parameter appears only in the combination J2<Pi ~ @- The 
minimum of the potential is given by the solution to the following equations: 

. ,9-TAi, N c Mi . , 

sm( — — ) = ___sm0 f . (5) 
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In the limit where all the quark masses are equal and Mj -C E, the approximate solution 
is given by 0j ~ In particular, in this limit and for Nf = 2, the vacuum energy V vac {6) 
can be approximated as follows: 



which was used in the Introduction for estimates. From this we see that a # irid -state is 
degenerate with the 9 = state in the chiral limit where m q = 0, as expected. 

3 Induced 6- Vacua 

It was recently argued [|]] that it may be possible to create an arbitrary 9 ^ state in 
heavy ion collisions. In this scenario, we propose that the # m<J -state is separated from our 
world (9 = 0) by a shell of debris expanding out at the speed of light. This idea is similar 
to the so called disoriented chiral condensate that has been extensively studied in the last 
ten years 

We would like to stress the point that an induced ^-parameter is very different from 
Qfund yjYAoh is zero in our world and which cannot be changed. The simplest way to 
visualize 9 md is to assume that right after the QCD phase transition the flavor singlet 
phase of the chiral condensate is non-zero in a macroscopically large domain. This phase 
is identified with 9 md . This identification is a direct consequence of the transformation 
properties of the fundamental QCD Lagrangian under U(1)a rotations by which the chiral 
singlet phase can be rotated away at the cost of introducing 9 md . From now on we will 
refer to a # md -state as simply a #-state, omitting the "ind" label. 

Numerical lattice calculations in the past have provided strong evidence that QCD 
undergoes a phase transition at a temperature in the range of 150—200 MeV. It is believed 
that above this critical point that chiral symmetry is restored and/or deconfinement occurs 
(quark-gluon plasma). The idea of the Disoriented Chiral Condensate (DCC) arises when 
we consider what happens immediately after the phase transition upon cooling within 
heavy ion collisions. DCC refers to regions of space (interior) where the chiral condensate 
points in a different direction from that of the ground state (exterior), and separated from 
the latter by a hot shell of debris. The difference in energy between the created state and 
the lowest energy state is proportional to the small parameter m q (see Eq.@) for both 
the DCC and the 6*-state and is therefore negligible at high temperature. 

First, we recall the DCC scenario as given in || with emphasis on the analogy between 
DCC and # md -state. Rajagopal and Wilczek [0] use the 0(4) linear sigma model to 
describe the low energy dynamics of the pions and the chiral condensate. All fields are 
represented by a 4- vector with components <fi = (o~, 7?), where a represents the chiral 
condensate and n represents the triplet of pions. Throughout the exterior region the 
vacuum expectation value of <fr is (v, 0). In the interior region, however, the pion fields 
can become non-zero and (cr, 7?) wanders in the four dimensional configuration space. 
The high energy products (shell of debris) of the collision expand outwards at relativistic 
speeds and separates the misaligned vacuum interior from the exterior region. 



(0\jtGG\0) 



V vac (9) ~ K 



vac 




(6) 
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In H all calculations were done under the assumption of a quenched system. In the 
quenched approximation, the (a, 7?) fields are suddenly removed from a heat bath (T > 
T c ~ 200MeV) and evolved according to the zero temperature equations of motion. This 
is done numerically by giving the fields a non-zero vacuum expectation value (0) 7^ and 
letting this field configuration evolve in time according to the zero temperature equations 
of motion. The hope is that regions of misaligned vacuum with an arbitrary isospin 
direction will be created. In || it was shown that if the cooling process is very rapid and 
the system is initially out of equilibrium, there is a temporary growth of long wavelength 
spatial modes of the pion field corresponding to domains where the chiral condensate is 
approximately correlated. In the case of DCC, the created state will relax to the true 
vacuum by coherent emission of pions with the same isospin orientation producing clusters 
of charged or neutral pions. To be more specific, let us consider the case Nf = 2. The 
matrix U is parameterized by the misalignment angle and the unit vector n in the 
isospin space: 

U = e **<M ,UU+ = 1, (¥ L ¥ R ) = -\(* L itf R )\ Uij. 
The energy density of the DCC is determined by the mass term: 

= -\Tr{MU + MW) = -2m|(M)|cos(0). (7) 

Eq. (0) implies that that any 7^ is not a stable vacuum state because ^ 7^ 0, i.e. the 
vacuum is misaligned. Since the energy difference between the misaligned state and the 
true vacuum is proportional to m q , the probability to create a state with an arbitrary 
at high temperature T ~ T c is proportional to exp[V(Eg — E )/T] and depends on only 
very weakly, i.e. is a quasi-flat direction. Right after the phase transition when fifty) 
becomes non-zero, the pion field oscillates. 

Now we turn to our main subject, the ^-vacuum state. First we would like to com- 
ment that even though the 9 parameter is related to the large constant E in Eq. (|J), 9 
dependence is actually proportional to m q since in the chiral limit all dependence can be 
removed by performing a U(1)a rotation of the chiral condensate phases. When the com- 
bination of (X}0i — 9) is close by an amount 0(m q ) to its vacuum value, the Boltzmann 
suppression due to the term E at high temperature is absent. This is essentially what 
allows the induced 6*-state to be formed. 

To model the creation of an arbitrary #-state, we will use the effective chiral Lagrangian 
described in the previous Sect. 2. We consider three flavors of massive quarks, and derive 
the following equation of motions: 

(pi + V 2 0i + E sin( ^ "L^ — ) - Mi sin(0j) - yfa = 0, i = u,d,s, (8) 

where the constant 7 in front of represents the damping of the system which can be 
attributed to various sources such as expansion or the emission of pions. 

In order to model the situation numerically, several assumptions are made. First, 
assume that the 9 md parameter acquires a non-zero value when the temperature T ^> m q 
and the system is out-of-equilibrium. Once again, 9 md 7^ follows from the assumption 
that the mean value of singlet phase (J2 4>i) is not zero over a macroscopically large domain 
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12 3 4 



Time (10"^ sec) 

Figure 1: \<f>i(k = 0)| is plotted as a function of time for the up, down, and strange quark. 
Notice that the zero momentum modes of the </>j fields settle to a non-zero value in a time 
on the order of 1CT 23 s. The times t\ and t 2 represent the values we chose for r s h e ii, the 
the time when the shell separating the two regions disappears. 



immediately upon cooling through the phase transition. Performing a U(1)a rotation then 
gives us an equivalent description with 6 md ^ 0. The second assumption is that the phases 
(j) u , <pd, and 4> s have small random values (modulo vr/3) at these high temperatures. The 
final assumption is that the cooling process is very rapid and the system obeys the zero 
temperature equations of motion after this initial quench. There is obviously some sort 
of damping mechanism present in any real system, and this damping can be attributed to 
the expansion of the region inside the shell and/or the radiation of Goldstone bosons. The 
damping constant, 7, was chosen to be of the same order order of magnitude as Aq C d, 
which should be a reasonable value according to all other scales in the system. 

The numerical algorithm used to solve these coupled differential equations was a 
fourth-order Adams-Bashforth-Moulton predictor-corrector method. The phases were ini- 
tially assigned random values according to a uniform distribution (|</>j| < 7r/3). The phases 
were rescaled in such a way as to be dimensionless. We evolved the equations for 6000 
time steps of size 10~ 27 sec on a cubic lattice with spacing a = 0.005 MeV^ 1 = 1 fm. 
In order to examine the momentum dependence, the Fourier transform [pjl of the field 
configuration was calculated at periodic intervals in time. The data was then binned 
according to the magnitude of the wave vector k so that the dependence on \k\ could be 
examined. We consider a volume of (8 fm) 3 . In all calculations that follow, we use the 
current quark masses m u = 4 MeV, irtd = 8 MeV, and m s = 150 MeV. 

In Fig. [I] we show that the zero modes of the 4>i fields settle down to a constant value 
in a time r ~ 10" 23 s. In the same time interval, all other modes die off. This, along with 
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a test of volume independence, shows that a true non-perturbative condensate has been 
formed. In the case where all quark masses are equal, the zero mode is expected to settle 
down to <f>i ~ O/Nf, which has also been verified numerically. 

In order to incorporate the effects of expansion, we could introduce a time dependent 
lattice spacing a(t). If a(t) increases with time, this will model the dilution of energy 
in the system. Following an analogy with an expanding universe, the coefficient of the 
damping term (f>i, 7 goes like constant x a /a. Although we have not taken this into account 
in our analysis, we can estimate the result as follows. The Laplacian term in Eq.(^) is 
numerically proportional to 1/a, and if a(t) increases with time, the system will require 
more time to relax to the 0-state. If the coefficient, 7, is altered as mentioned above and 
a(t) increases at a constant rate, the damping term will become smaller as time goes on. 
The overall effect of a time dependent lattice spacing is to increase the formation time of 
the 0-state. 



4 Goldstone Bosons 

Now that we have argued the induced #-state may be formed in heavy ion collisions we 
would like to discuss how to detect them experimentally. Naturally, the first thing to 
look for would be axion production. If axions[|l4| do exist, they can be produced during 



the relaxation of 8 md 7^ to the lowest energy state with 8 = 0. However, the axion 
production rate under the conditions which can be achieved at RHIC will be too low 
in comparison with the limit already achieved from the astrophysical and cosmological 



consider at ions JT^]. Therefore, the question remains: "How could the the formation of a 
#-state be detected experimentally?" Fortunately there are a number of ways in which 
this might be done. One proposed possible signature is related to the observation of 
specific CP odd correlations of the Goldstone bosons but the effect would probably 
be washed out by final state interactions. A second possible signature was proposed in 
II which comes about because of the unusual properties of the Goldstone bosons and 
rf in the ^-background. In the #-world the masses of these particles are altered and for 
example the mass of 7r° would be noticeably smaller. As well, these particles become a 
mixture of pseudoscalar /scalar rather than pure pseudoscalars as for 8 = 0. This fact 
means that the CP odd decays n — > nn and rf — > mi are no longer suppressed in the 8 
background and become of order 1 f|. The full width of the n', n decays in our world are 
much smaller than the above decays in the 6*-world: 

r^otV -> vttt) ~ 2 MeV > T^(rj') -0.2 MeV (experimental), (9) 
Te^oiv -> tttt) ~ 0.2 MeV > T^oiv) ~ H8 keV (experimental). (10) 

Thus observation of an increased rate of decay of n', n at masses slightly shifted from the 
accepted values would be a signature of the ^-vacuum. However, this signal may be a 
difficult one to distinguish in the large amounts of data generated in heavy ion collisions. 
Therefore we would like to discuss a signature that would be easier to experimentally 
observe. In this work we consider an additional new signature which could possibly verify 
if a #-state is created in heavy ion collisions. As was discussed briefly in the introduction, 
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the creation of a 0-state could greatly enhance the production of low momentum (~ 
10 MeV - 100 MeV) Goldstone bosons. 

Before going into the details we would like to discuss the concepts behind this idea. 
We have already presented some arguments that a #-state may be formed in heavy ion 
collisions if protected from the exterior 9 = world by the shell of debris. At some point, 
however, this shell ceases to isolate the interior and the influence of the exterior world 
will be felt. The way in which this happens is not well understood, but we would like to 
use an instantaneous approximation in order to obtain some feeling about the possible 
outcome. There are basically two time scales present in the system: the lifetime r s h e ii 
of the shell separating the two worlds and the scale associated with strongly interacting 
particles, (Aqct>) _1 . If these differ by orders of magnitude, the disappearance of the shell 
can be considered either as an instantaneous perturbation or an adiabatic perturbation. 

Realistically, we expect T s h e ii to be about the same order of magnitude as (Aqcd)" 1 and 
therefore the disappearance of the shell is somewhere between an instantaneous process 
and an adiabatic one. However, we do not yet know how to treat the process properly 
so we will use the instantaneous approximation in order to obtain an order of magnitude 
estimate of the spectrum of emitted particles. 

The instantaneous perturbation is an approximation that is well known in quantum 
mechanics. For our specific problem we essentially assume that all states which were 
formed in the 9 %nd background will suddenly find themselves in the new vacuum state 
(with 9 = 0) after the shell instantaneously blows apart. These states have no choice but 
to transform to the asymptotic states of this new (for them) 9 = world. The dynamics 
of this strong transformation is quite complicated and beyond the scope of the present 
work. However, if this transition is instantaneous, the procedure to do the calculations 
for such a transition is well known: we assume that the initial state has no time to adjust 
its properties when the shell decays, and therefore, we expand the initial state in terms 
of the asymptotic states of 9 = world. 

In more detail, we do the following. The field values, (pi, obtained for the #-state as in 
previous work [||] are embedded into a larger grid where the field values take their 9 = 
vacuum values. This data now contains a plateau of field values corresponding to #-state 
values surrounded by zero field values. This field configuration must now resolve itself 
into asymptotic free particles. With this process in mind we determine the momentum 
spectrum of free particles corresponding to this distribution by considering the Goldstone 
boson fields as true quantum fields at this instant. 

In order to calculate the spectrum, we must obtain the quantities a(k) and a^(k). 
These are obtained by performing the Fourier transform of the field configuration at time 
t = r s h e u when the shell breaks down, 

a(k) = J d 3 xip(x) exp(— ik ■ x), (11) 

where ip(x) is the distribution amplitude of the <pi fields obtained from evolving the equa- 
tions of motion (Eq.(|S|)) on a cubic lattice, as was done in [|TJ. From Eq.(|IT]), the number 
operator is given by N(k) = Noa(k)a^(k) where Nq is an overall constant that will be 
determined using the following argument. As we mentioned earlier, the #-state differs 
from the true vacuum state by a small amount of energy ~ m q . The amount of energy 
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that is available when the #-vacua decays is obtained by analyzing the 9 dependence of 
the vacuum energy density E: 



E e = m q \(^)\N f cos(^-). (12) 

Therefore, the amount of energy AE V available due to the formation of the 0-state as 
discussed above is given by 

AEV = (E 0=o -E e ).Vc20 (p^s) GeV, (13) 

where V is the the volume of the created ^-region. If we take for example a #-state with 
volume V = (8fm) 3 the amount of energy available is AE V = 10 GeV. This represents 
only a small amount of energy compared to the cm energy, y/s = 40 TeV, expected in 
heavy ion collisions. If the total amount of available energy is known, the constant N 
can be fixed by enforcing the following conservation of energy constraint 

/d k — * j. — * 
^y 3 N oai (k)4(k), (14) 

where the sum is performed over all types of particles considered (ie. all neutral Goldstone 
bosons). Once this constant N is fixed, the total number of pions and etas can be 
calculated by 

N = I^k-^)^), (15) 

for the simplified case of only one collision of two nuclei. If we consider that RHIC will 
collide gold nuclei at a rate of central collisions of about 1 kHz, we can estimate the total 
number of these low momentum particles produced by simple multiplication. In what 
follows we present all results for just one event. 

As we have already described, we proceed as follows. At some time t before the fields 
4>i settle down to the constant field configuration (see Fig. [I]), we take the position space 
data and embed this in a larger square grid where the field values are zero (our world). 
We assume an instantaneous perturbation and then take the fourier transform of this 
field configuration in order to obtain the operators a{k) and a'(k). It may be argued 
that the instantaneous approximation is not completely justified in this case as we expect 
i 'shell ~ tqcd and therefore this calculation gives us a rough estimate at best. 

In order to analyze the spectrum of the diagonal (neutral) components of the matrix 
U (7T°, r), r)') we must analyze the following combination of the 

V = A^u + fa-Ws), (16) 
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For all numerical calculations, we can take the constant /,,/ to be ~ /„.. One of the crucial 
decisions which must be made is when to choose the time at which the shell breaks down. 
We chose several values for r s heii- The majority of the particles which are produced due to 
the formation of the |#)-state have momentum k < 25MeV ( for our choice of the volume 
V = (8fm) 3 ). In Fig. |2| we show a plot of N n o(k) as a function of \k\ for the neutral pion. 
The function N 7T o(k) has dimensions MeV , so that the total number of particles, A^o, 
is a dimensionless number given by: 

/oo 
dkN T o(k) 
-oo 

00 47T — * 1 — * 

dk N a(k)a\k). (17) 

-00 2,{2,TryWk 

The two different lines represent the times at which we assumed the shell to break down. 
The solid line represents the earlier time {r s heii — ti — 1500/6000 time steps) while 
the dotted line represents T s h e u = t<i = 3000/6000 time steps (this representation for t\ 
and t2 will be followed in all graphs that follow). For this calculation, we assumed the 
created #-state had a volume of (8 fm) 3 which exists in a larger volume of (64 fm) 3 . The 
explanation of this phenomenon is simple: at time t± the amplification of the zero mode 
is not as profound as at time ti. Furthermore, at time t\ a considerable portion of the 
energy goes to high momentum particles. 

Below we show the same graph for the 77 and the 77'. Notice that that the r/ is produced 
in large amounts compared to the 7r°-mesons and the r/-mesons. In order to obtain the 
total number of each particle produced, we use Eq.flTB"!) to calculate the total number 
of particles produced per collision: A^o = 0.1, N v = 2.8, N n i = 19.3 for t = ti and 
N n o = 0.5, N v = 4.3, Nrf = 18.4 for t = £2- We have also checked that in the limit 
where all quark masses are equal, only the 7/ is produced, as expected. It should be 
noted that, as demonstrated by these figures and the coarsening phenomenon (i.e. the 
phenomenon of amplification of the zero mode as time increases), we would expect that 
as Tgheii increases up to a maximum value, the majority of the particles would reside in 
the low momentum regime. We consider our estimate as a very conservative one due to 
the fact that we do not count the particles 77 and 77' which will be produced during the 
formation period of the 9 tnd state. These particles will have higher momentum modes 
and can not be easily distinguished from the large number of the similar particles from 
the background. Therefore, the main conclusion of our analysis is as follows: if induced 9 
state is produced than it will result in the strong enhancement of emission of low-energy 
77 and 7]' mesons which clearly signals the existence of the induced 9 state. 

When our analysis was completed, we became aware that a similar conclusion has 
been reached in a recent paper by Baier et al. [PJ where the production of an excess of 
77'-mesons is discussed in the scenario suggested in J7| . Despite the fact that the starting 
points of the present work and [T(| are quite different (we start from the induced 9 



state while in paper [TIJ the starting point is a metastable state which acts like a region 
with a non- vanishing 9 angle), the general conclusion is very similar. Namely, a strong 
enhancement of emission of low-energy T/'-mesons is expected irrespectively of the origin 
of the initial unusual state. The details, however, are quite different. The main difference 
is the following. We assume that a 0-state is produced and instantaneously decays such 
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Figure 2: The number of neutral pions produced, N n o(k), is plotted as a function of the 
magnitude of the wave vector, \k\. Above we show that the momentum distribution of 
the 7r°-mesons produced is primarily < 25 MeV for two different values of r s h e ii. The 
solid line represents the earlier time t\ and the dotted line represents the later time £2 (see 
Fig. [I] for the positions of t\ and £2 in the evolution of the (pi fields) . This is expected as 
the formation of the #-state can be attributed to the enhancement of the low momentum 
modes. 
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Figure 3: N V (K) is plotted as a function \k\, at two different times t\ and as in Fig. 
|2|. The higher peak represents the later time with the main difference being that now a 
larger percentage of the produced particles lie in the A; < 25 MeV range. 
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that we calculate the particle production due to this single event, not taking into account 
previous particle production during the formation period. As was mentioned above, we 
expect that those particles will have much higher momentum modes and we ignore them. 
The approach of |16| is different from ours as their calculation of the number of produced 
7/-mesons is based on the specific mechanism of the decay of the original metastable state 
0. Such a decay receives contributions from the whole evolution, not just one particular 
instant. This difference explains the difference in rates: for a domain with radius 5 fm, 
they estimate that about 90-100 r/-mesons would be produced while our estimation is 
4 — 5 times smaller. 

In the above calculations of N(k) the geometry of the created 9- state is assumed to 
be a cubical box. For the case of heavy ion collisions, the more realistic case is that of an 
ellipsoid or elongated rectangular box. Upon collision, we expect the collision to create 
a region of quark-gluon plasma which has an asymmetry in one direction. In order to 
model this situation, we evolve a rectangular grid and embed it in a larger square grid 
and compute the Fourier transform. Once again, we consider the angular averaged value 
ofiV(fc). 

For a ^-region with dimensions 32 fm x (8 fm) 2 , we show the spectrum of the rf at 
the same times as shown in Fig. The side of length 32 fm is the direction parallel 
to the beam direction. The low momentum modes still dominate, but the peak is not as 
sharp and the higher modes show a stronger presence. Also, since the ^-region is now 
larger, from Eq.flTB"!) we would expect more particles to be created since there is now more 
energy available. This is evident when comparing Fig. [| with Fig. [|. For the two different 
times shown, we find that N n o = 0.3, N v = 12.7, N v > = 76.3 for t = t\ and N n o = 1.7, 
N„ = 16.8, Nrf = 73.8 for t = t 2 . 
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Figure 5: The momentum distribution (N v '(k)vs.\k\) is now shown for the case where 
the geometry of the ^-region is considered as a rectangle. We consider the case where the 
size of the ^-region is larger along the longitudinal directions compared to the transverse 
directions. The solid and dotted lines once again represent times T s h e ii = h and t% respec- 
tively. Although the majority of the spectrum still exists in the low momentum range, 
the higher momentum modes now become more evident. 



In this case where spherical symmetry is lost, we must also examine the angular 
distribution of the system. We choose the coordinates such that the long edge of the 
^-region is parallel to the z-axis. In Fig. |6|, we show the angular distribution of N n o(6) as 
a function of the azimuthal angle 9. This plot suggests that the majority of the particles 
would have larger momentum components in the direction perpendicular to the beam 
direction. 

As was shown in Figs. HHlJj], and ||, if a #-state can be created in heavy ion collisions we 
should expect an excess of low momentum Goldstone bosons. In particular, the r/ would 
be produced in large amounts. Besides that, the low momentum 7r°, rj, and //-mesons 
would considerably increase the photon and e + e~ pair production through such decays as 
7t° — > 7e + e~, rj —>■ 7e + e~, and rj' — > / ye + e~ . In particular, it would result in enhancement 
of low-energy dileptons, which could possibly provide a solution to the observation of an 
excess of these particles seen at CERN [[J. Indeed, as was demonstrated by Kapusta et al. 
in (T7| a sufficiently large enhancement^ of ^'-production could easily explain the excess 
of dileptons seen at CERN || . The large amount of produced 7/-mesons in our work is a 
direct consequence of the decay of the induced #-state. 

2 The enhancement mechanism of //-production suggested in |l7j] is quite different from what we have 
considered in this paper. However, the general conclusion (that rj' enhancement leads to the excess of 
dileptons) is insensitive to the nature of the 7/-mesons. Only by performing a detailed experimental 
analysis of the momentum and angular distributions can one mechanism be chosen over another. 
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Figure 6: For the scenario shown in Fig. [5], we consider the angular dependence on the 
azimuthal angle for the 7r° instead of the 77'. The azimuthal angle is defined so that 6 = 
coincides with the beam axis. 

5 Conclusion 

In this work we showed that the consequences of creating non-trivial #-vacua in heavy 
ion collisions could be the amplification of production of light Goldstone bosons in the 
10 MeV momentum range. In all calculations, we worked in the instantaneous approx- 
imation where the shell separating the (9-state disappears in a time much less than any 
internal time scale. These low momentum particles could possibly then decay to low 
momentum photons and dileptons, which could be easily detected. We would like to 
make the suggestion that this could possibly account for the unexplained abundance of 
low momentum dileptons observed at CERN M. We are thankful to Dima Kharzeev for 
correspondence regarding the relevant work in [ |17|j . 
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